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Abstract 



i, 

t^ , A new theoretical analysis of the current-charge and charge-charge propaga- 



tors is presented for the Hubbard model, using the static approximation for 



'^ I the Composite Operator Method. This approximation manifestly preserves 

O I a sum rule which governs the single-site dynamics. We compare our results 



with those obtained by numerical analysis. 
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1. The study of optical spectra is an important probe of the dynamics of strongly 
correlated electronic systems and gives information about the normal and superconducting 
states of these materials 0,0] • Also, the optical conductivity is of fundamental theoretical 
interest because the spectral weight at low frequencies might be a natural order parameter 
for the Mott transition ^. 

The optical spectrum of undoped cuprate compounds, such as La2CuOi, Y Ba2Cu^0Q 
and Nd2Cu04, is characterized by the 02p — Cu3d^2_y2 charge-transfer excitation with a 
threshold energy in the 1.5 — 2.0eV range |l|,§]. By doping charge carriers, a drastic change 
in the optical spectrum occurs. The charge-transfer excitation gap is rapidly suppressed 
and instead an edge shows up around l.Oey in the reflectivity spectrum, indicating that 
the spectrum is dominated by low-energy excitations in the energy range lower than l.Oel^ 
0. This kind of changes with carrier doping are universal for all known hole-doped and 
electron-doped cuprate superconductors. Besides, the appearance of optical anomalies is 
not unique to the doped copper oxides, but is almost universally observed in the spectra of 
the systems with strong electron correlation ljl],0,|§|. 

In the simplest form, the Hubbard model is believed to describe many properties of 
strongly correlated fermion systems. The applicability of the model to the superconducting 
copper-oxides is related to the fact that lower and upper Hubbard subbands can mimic the 
02p and CuSd^^^yi bands of the cuprates and thus the charge-transfer gap can be described 
by an effective Hubbard on-site coupling constant. 

Theoretical studies of the optical conductivity for the Hubbard model started long time 
ago : it was investigated by the moment method, by the equation of motion method for 
the Green's functions in the Hubbard I approximation and in the Hubbard HI approxima- 
tion. For the two-dimensional Hubbard model the strong coupling limit has been considered 
mostly by numerical methods based on exact diagonalization for small clusters P-p^. A 



detailed discussion of optical and photoemission sum rules for one- and two- dimensional 
Hubbard model has been given recently by comparing a perturbation theory in the kinetic 
term with numerical calculations |l^ . 



An analytical analysis of the frequency dependent conductivity has been given by per- 
turbation method in weak-coupling limit |T6| , p^ and by applying the memory function tech- 
nique in terms of the Hubbard operators ||^. Recently the optical conductivity has been 
investigated within the djTiamical mean- field approach [6] . However, for the realistic planar 
Hubbard model non-local corrections to the transport vertices and the self-energy are impor- 
tant [6]. Beyond this, vertex corrections are crucial to satisfy the relevant Ward-Takahashi 
identities [17]. 

In this Letter, we present a new theoretical analysis of the conductivity tensor by use of 
the Composite Operator Method [18] for the two-dimensional Hubbard model. By exploiting 
sum rules after the Ward-Takahashi identities and the Pauli principle, the calculation of the 
current-current propagator is reduced to the charge-charge one in the framework of a fully 
self-consistent scheme. Further, it will be shown that the charge-charge propagator preserves 
an important sum rule which governs the single-site dynamics avoiding double occupancy 
configurations unless to form local singlets. 

The Letter is organized as follows. First, the model is set up and the formula for the 
conductivity tensor in the linear response theory is summarized. The Ward-Takahashi iden- 
tities are given in the context of the two-dimensional Hubbard model and the hierarchy 
degrading the current-current propagator to the current- charge and charge-charge ones is 
built up. The calculation of the current-charge and charge-charge propagators is realized 
within the static approximation for the Composite Operator Method. Then, the results 
for the optical conductivity and the Drude weight are presented and compared with those 
obtained by numerical analysis of small clusters. Some concluding remarks are presented at 
the end. 

2. In the presence of an electromagnetic field described by the vector potential A{r,t), 
the Hubbard model is as follows 

H = Y: lUM) - ^^5^,] ^H^J) + UJ2n^i^)n^i^) (1) 

where we use a standard notation. The hopping matrix elements contain the Peierls factors 



which guarantee the gauge invariance. 

In the framework of the hnear response theory, the electric conductivity tensor is given 
as follows 

(TabiK^) = e^- — -—-{Ka)6ab " T7 — -—-gabCk, Uj) (2) 

t[uj + tr]) i[uj + 17]) 

where Qahi}^,^) is the retarded current-current propagator, with Kaii) being the kinetic 
energy density operator in the a direction. 

In particular, the frequency-dependent uniform, i.e. k = 0, electric conductivity a^xi!^) 
is given by 

ai{uj) = Rea^xi^) = D5{u) + Oindi^) (3) 

where D is the Drude weight 

(4) 



D = e^7r 



-(ir,) + -Re^,,(0,cj^O) 



e2 



which measures the ratio of the density of the mobile charge carriers to their mass [3] , and 
(Tine is the incoherent contribution 

o"inc(^) = Im5f^^(0,t^) (5) 

UJ 

By defining the charge p{i) and current j(i) densities as 

(6) 
c{^) (7) 



p{i) = ec^i)c{i) 






-V 



it is immediate to obtain the conservation law 

V-j(^) + ^p(0 = (8) 

The symmetry content of the algebraic equation (H) manifests at level of observation as 
relations among matrix elements once a choice of the physical space of states has been made. 
Indeed, by defining the causal charge and current propagators as Xabihj) = {T [ga{i)gb{j)]) 



where ga{i) = {pij), jxif), jyij)) for a = 0, x, y, we can derive a series of Ward-Takahashi 
identities, which in momentum space read as 
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(9) 
(10) 

(11) 

(12) 
(13) 



By these identities we see that the current-current propagators Xa6(k, cj) can be expressed 
in terms of the current-charge and the charge-charge ones. 

In order to calculate the current- charge and the charge-charge propagators, we use the 
static approximation [19] for the composite field 



m 






(14) 



By defining 



N{k) = F.T.{T[n{{}n{j)]) S{k) = F.T.{T[^{{}^^{j)]) 



F.T.{{^(t),^\j)}^^^) 



(l 



- 

2 " 



2/ 



(15) 



(16) 



we have 



N{k) = -in'{2Tifa~H^''\k) 



^^3^{/n [Qiiii(A;) + giii2(A;) + gi2ii(A;) + gi2i2(A;)] 

+ /22' [Ql212{k) + Ql222{k) + Q2212{k) + Q2222(A;)]} (17) 



with 



iO, f 



(18) 



where D is the double occupancy 



D = ^{v\i)v{^)) 



(19) 



k is the four-dimensional vector k = (k, cu); the symbol F.T. denotes the Fourier transform. 
Then, it is direct to see that 



Xoo{k) = e'N{k) 



(20) 



and 



Xaoik) = (-zte2a)(e^^- - l)N{k) {l^,' [(C^^) + {Cv^)] - I,,' [(r/"e^) + (r^"r^t)] } 

-2I^^\-tte^a) [Ql,,,{k) + Q1,^,{k) + Ql^,^{k) + Ql,^,{k)] (21) 



where we have defined 



QIb,6%^) 



iVL 



cfpdn 



./37.5 V"' --/ (^27r)3 

We have in the coordinate space 



pi{pa+ka)a _ p-ipaa 



S^(s{K^)S^5iP + K^ + ^) (22) 



N{i,j) =n 



n{2 — n) 



{in iQuiiihj) + Quuihj) + Quuihj) + Quuihj)] 



n-2D 

+-'22 [Wl212v^'.^J ~r Wl222v^'iJ "•" W2212\^'iJ "•" n;2222 V^' i Jj J 



In particular, at equal sites (53) reduces to 



N{i, i) = {n{i)n{i)) = n + 2D 



(23) 



(24) 



We thus see that the charge propagator satisfies the sum rule 



(2^ 



/ d^kduxooik) =e^{n + 2D) 



(25) 



as it should be. The physical content of (|25|) is that of the Pauli principle because it governs 
the single-site dynamics of two fermions. 

By means of the Ward-Takahashi identity (p!2| ) we obtain 

2 [1 — cos(Ka;a)J 2 

where K = — 4t(c^(2)c"(i)) is the kinetic energy per site. From the analytical structure of 
equation (p6|), we note that Xxx{kx,ky = 0,0) = o?e^K/2 does not depend on k^ and is 
proportional to the kinetic energy. It is remarkable that such a peculiar behaviour has been 
previously found by use of Quantum Monte Carlo techniques [8]. The correspondence with 
the notation of Ref. 8 is given by g^x = ~^^xx- 

By taking the limit A;^. ^ in (p6D and by considering the retarded function we have 

,..(0, u;) = -a e ir + hrn^ , ^, _ ^^^^.^^^j (27) 

Recalling (|^) we have for the optical conductivity 

ai(cu) =L'(5(cu) + (Ti„,(cu) (28) 

where the Drude weight D is given by 

D = vrRelimlim4^^^!^i%4 (29) 

^^Ofc,^o2[l -cos(A;^.a)] ^ ' 

and the incoherent part OinciyS) is as follows 

OincK^) = Im hm — -—— (30) 

uj k^-^o 2[1 — cos(A;2;a)J 

As emphasized in the beginning, the whole problem of the response to an external electro- 
magnetic field has been reduced to the evaluation of the charge-charge propagator. In the 
static approximation of the Composite Operator Method [18-20] the charge-charge propa- 
gator, as the spin-spin one, can be connected to convolutions of single-particle propagators 
[19]. It is important to note that the convolutions in ( PT] ) and ( P^ involve higher order 



single-particle Green's functions. Then, such a scheme does not contradict the conclusion 



reached by exact diagonalization analysis [11] that any attempt to compute the conductivity 
from a convolution of single-particle propagators is destined to fail in one dimension and is 
questionable in planar systems. In our case, the occurrence of convolutions is related to a 
linearized dynamics, but also involves choice of occupation dependent electronic excitations 
as basic fields [19]. In this way, the electron propagation is described as a repetition of 
composite excitations which automatically take into account scattering of electrons on spin 
and charge fiuctuations due to strong correlations. Also, the use of a higher order basic field 
gives the advantage of implementing the Pauli principle and of a clear theoretical under- 
standing of the terms originating from intraband or interband propagation. In fact, it turns 
out that only intraband excitations contribute to the Drude weight, whereas interband ones 
build up the incoherent part. 

In the static approximation, by means of the previous results, the retarded charge-charge 
propagator is given by 

n{2 — n)e^ 



xW^,oo) = -^n'(^\2nfa-H^'\k)- 



n-2D 



{in [Qnuik) + Q^U2ik) 



Ql21l{k) ~^ Quuik) +-'22 Ql212(^) + ^12221^) + ^22121^) + ^22221^) J \'^^) 



where 



Q%^siKu^)= E 



n 



«j= 



i(2 



n 



<iV^i^!^i4^T^i^^.S(p)-S'(k + p) (32) 



Lj + Ei{p) -Ej{k + p) +ir] 



f{uj) is the Fermi distribution function. cTaUp) are the spectral functions of the single- 
particle propagator. These functions have been previously calculated [18, 19]. 
Let us define 



A,(k, p) = In' [cri?(p)4^(k + p) + a«(p)aii)(k + p) + ag(p)aii^(k + p) + a!?(p)ai?(k + p) 

ag ip)ag^ (k + p) + aS (p)4^ (k + p) + 4? (P)^i'2^ (k + p) + 4? (P)4? (k + p)! (33) 



-'22 



and 



Xij{k,uj) 



Q 



(2vr) 



d^ 



f[E^{k + p)]-f[E,ip)] 



UJ 



E^ip) 



E,(k 



P) 



11] 



-Aij(k,p) 



(34) 



Then, x^(k, a;) can be rewritten as 
and 



D = _<^-^>^' V lim lim ^^^^R^-^-#->^) (36) 

It is direct to see that the interband terms Xi2(k, cj) and X2i(k, cj) do not contribute. Also, 
by means of straightforward calculations the contribution of the intraband terms gives 

""Ji'-.TTTTnr = -WtWv f^P-->P^)UP) (37) 

uj^Ok^^o 2 [1 — cos(A;2:a)J 4A;bJ (27r)"^ j 

and the Drude weight takes the expression 



where we defined 



i(o) 



U^)=<{v)[l-mp)W{v) (39) 



with 



,,p,.,,„,(^M^) (40) 

In the same way, it is straightforward to show that the incoherent part is given by 

ln(2-n)e2 2 a'^uj'^Xij{k^,uj) 

ainc{uj) = — — 2^ Im hm — '-— — - (41) 

u n-2D .^^ k^^o 2 [1 - cos(A;:ra)] 

where the intraband terms do not contribute. In particular, it is found that o"j„c(^) = 
for uj < uJc- For any finite value of U there is a gap in the incoherent part of the optical 
conductivity, given by 



Uc = 2^/n/22f/ [At half-filling uJc = U = 2fi] (42) 

In Fig. 1 we present results for the Drude weight and from numerical analysis of small 
clusters [13]. 

9 



For f/ = 15 the Hubbard gap is clearly visible at half-filling because, as discussed by 
Kohn [3], one expects that in an insulating phase D will vanish. Also, there is no spectral 
weight in the interval < cj < 15 [cfr. Eq. (^21)]. Upon doping, spectral weight is transferred 
from the region above the Hubbard gap to the Drude peak at zero frequency. This is up to 
a critical doping where there is a deviation downwards which may be taken as a signature of 
a drastic change in the nature of the carriers. Indeed, after this doping the system starts to 
resemble a gas of non-interacting electrons and the Drude weight follows the kinetic energy 
[12]. It is worth noting that the value of the Drude weight and of its slope in proximity of 
half-filling should be carefully interpreted within the scheme used. In fact, many numerical 
and theoretical indications [2] signal the onset of additional states in the gap due to doping, 
extending upwards from the lower Hubbard band, which absorb a large spectral weight. 
Then, the absence of a mid-gap band accompanied by a conserved total spectral weight, as 
in the approximation used [20], will very much affect the spectral weight transfer between 
the interband charge excitations and the intraband ones. Further, this occurrence can be 
responsible of a Drude weight larger then the non-interacting value near half-filling for some 
particular choice of the on-site Coulomb repulsion U . 

The kinetic energy is shown in Fig. 2. The results have been obtained through the 
scheme given in Refs. 21 and 22 which makes extensive use of thermodynamic identities for 
the Helmholtz free energy per site. 

These results to a large extent contradict data available both from numerical and exper- 
imental analysis [1, 2]. The observed ratio uj^/U , connected to the energy scale associated 
with interband charge excitations, is greater compared to what is expected from Lanczos's 
techniques [12]. uOc/U is depicted as function of the filling n in Fig. 3. 

However, the small system sizes and the lack of theoretical analysis of finite size cor- 
rections make interpretation of the results uncertain. A more troublesome discrepancy is 
related to the absence of a clear feature centered at low frequency beyond the Drude peak, 
quickly developing by varying the dopant concentration, as already discussed in relation with 
the Drude weight results. The presence of a mid-gap band is a genuine two-dimensional ef- 

10 



feet [11] that eannot be related to interband charge excitations which lead to an absorption 
band centered near U , as seen also in the strong coupling counterpart of the Hubbard 
model [2] (i.e., the t — J model). This is the well known mid-infrared band, which has been 
observed in several cuprate superconductors and strongly correlated electronic systems by 
use of analytical and numerical techniques [1, 2]. We do not at present have a resolution of 
this discrepancy, but a possible one is as follows. 

The energy scale on which the center of gravity of the mid-infrared band is centered 
is 4t^/f/, that is the value of the Heisenberg exchange J when derived from the one-band 
Hubbard model. Then, optical mid-gap effects can be observed only if the scheme presents in 
addition to the lower Hubbard band a quasi-particle peak separated by a gap or pseudogap 
of order J. On the contrary, we cannot reproduce the latter result since for our estimations of 
the single-particle propagator we have used a two-pole expansion for the one-electron spectral 
density. In this respect, an enlarged composite basis with a third asymptotic excitation can 
resolve this structure. 

In conclusion, we have shown that the operatorial equation (^ expressing the charge con- 
servation law has a direct corollary in terms of the existence of Ward-Takahashi identities, 
which hierarchically degrade the evaluation of the current-current propagator to the current- 
charge and charge-charge ones. Such an occurrence takes place also in the non interacting 
system (i.e., U = 0). However, it is worth pointing out the different role played by the Ward- 
Takahashi identities in non interacting systems compared to that in the framework of an 
approximate theory for interacting ones. In the first case the Ward-Takahashi identities are 
automatically satisfied and it is a matter of convenience to compute the conductivity from 
the current-charge or charge-charge propagators. In a theoretical scheme for an interacting 
system there is no way to avoid some approximation. Then, in addition to a clear conve- 
nience in computing a hierarchically degraded propagator, there is the possibility to build 
up a physical space of states where the general symmetry principle of charge conservation 
law can be enforced to model the dynamics. The scheme turns out to be even more coher- 
ent by preserving the fermionic character of the propagating charges or, in other words, by 

11 



implementing the Pauli principle. That is, by recovering such general symmetries the prop- 
agation of the charge and their intrinsic dynamics is constrained in a suitable Hilbert space 
where the charges propagate without degrading in number and avoiding double occupancy 
configurations unless to form local singlets. 

The authors would like to express their thanks to Adolfo Avella for valuable discussions 
and Sarma Kancharla for a careful reading of the manuscript. One of us (D.V.) thanks 
Gabriel Kotliar for kind hospitality during his stay at the Rutgers University where a part 
of this work has been done. 
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FIG. 1. The normalized Drude weight D„ = D/(27re^) as a function of the fihing n for various 
couplings U . Full squares are zero temperature data using exact diagonalization techniques on 
4x4 sites. Open squares and triangles indicate results for a \/lO x a/IO site cluster. 
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FIG. 2. Kinetic energy per site —K as a function of the filling n for various couplings U . Full 
rhombi are zero temperature data using exact diagonalization techniques on 4 x 4 sites. Open 
rhombi indicate results for a vlO x vlO site cluster. 
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FIG. 3. iOc/U as function of the filling n. 
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